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We study characters of recent type la supernova (SNIa) data using evolving dark energy models 
with changing equation of state parameter w. We consider sudden-jump approximation of w for some 
chosen redshift spans with double transitions, and constrain these models based on Markov Chain 
Monte Carlo (MCMC) method using the SNIa data (Constitution, Union, Union2) together with 
baryon acoustic oscillation A parameter and cosmic microwave background shift parameter in a flat 
background. In the double-transition model the Constitution data shows deviation outside la from 
ACDM model at low (z < 0.2) and middle (0.2 < z < 0.4) redshift bins whereas no such deviations 
are noticeable in the Union and Union2 data. By analyzing the Union members in the Constitution 
set, however, we show that the same difference is actually due to different calibration of the same 
Union sample in the Constitution set, and is not due to new data added in the Constitution set. 
All detected deviations are within 2a from the ACDM world model. From the ACDM mock data 
analysis, we quantify biases in the dark energy equation of state parameters induced by insufficient 
data with inhomogeneous distribution of data points in the redshift space and distance modulus 
errors. We demonstrate that location of peak in the distribution of arithmetic means (computed 
from the MCMC chain for each mock data) behaves as an unbiased estimator for the average bias, 
which is valid even for non-symmetric likelihood distributions. 

PACS numbers: 98.80.-k, 95.36,+x 



I. INTRODUCTION 

The type la supernova (SNIa) observation indicates 
current accelerated expansion of the universe [H,[2|. In or- 
der to explain the acceleration the concept of dark energy 
was introduced with many different realizations using a 
fluid, a field, etc. The SNIa data arc important to study 
possible evolution of dark energy, often parametrized by 
an equation of state parameter w = p/ pc 2 ; p is the pres- 
sure and pc 2 is the energy density of the dark energy com- 
ponent. Although the SNIa data have been continuously 
improving ; it is still remarkable that all current ob- 
servational data are well accommodated by the simple 
ACDM model with w = — 1. Thus, at the moment the 
observational data do not seriously demand evolution in 
the dark energy equation of state. 

Although ACDM model explains observational data 
well, diverse time-varying dark energy models and anal- 
ysis methods have been introduced to deal with increas- 
ing amount of data 
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Possible evolution of dark 
energy has been probed mainly by Chcvallier-Polarski- 
Linder (CPL) parametrization w(a) = WQ + w a (l — a/ao) 
[iol . |4~0 | and the redshift-binned parametrization method 
based on the decorrelating technique [HI, H3, HH . In the 
latter, the equation of state parameter in each redshift 
bin is estimated independently of others based the piece- 
wise constant, linear or cubic spline parametrizations. 
Compared with the CPL parametrization, which has only 
two free parameters with special evolution, the redshift- 
binning method is more suitable to know specific dark 
energy variations because we can divide the redshift in- 
terval into many bins as we want. Although there is a re- 
striction on the number of redshift bins due to insufficient 



data, the rcdshift-binning method is more appropriate to 
study the nature of dark energy. 

There are several works showing tendency of w deviat- 
ing from —1 especially in the low redshift region (32| - [38j . 
The tendency indicates slow down of acceleration with 
the deceleration parameter q increased in z < 0.3 [32| . 
Similarly, authors of Ref. [33| showed that w is rapidly 
decreased as z increases under both the CPL and the 
redshift-binned parametrization methods. In Ref. [HI, 
the linear spline method was used to show that w(z) can 
be oscillating around w = — 1 with two turning points 
where w < —1 at z ~ 0.4 and w > —1 at z ~ 1. The au- 
thors of Ref. [3{| also reported weak hint of possible evo- 
lution of dark energy; it is found that under the piecewise 
constant parametrization w is less than — 1 in z between 
0.25 and 0.75. On the other hand, no significant evidence 
for dark energy evolution is found based on the similar 
methods in Refs. [36l - [38j . In the high redshift region, 
authors of Ref. [3!| report that the dark energy deviates 
from the cosmological constant at redshift z > 0.5 based 
on analysis of SNIa and gamma-ray bursts (see also Ref. 

In this paper we study the evolution of dark energy us- 
ing a double sudden-jump transition model with varying 
transition redshift, which is a simple version of the piece- 
wise constant parametrization method. To constrain our 
model we use the recent SNIa data such as Union 0], 
Constitution 0], Union2 0] samples, together with the 
baryon acoustic oscillation (BAO) parameter [12] and 
the cosmic microwave background (CMB) shift param- 
eter [43l - [45l ] . Especially we investigate the characters of 
the recent SNIa data sets in view of the simple evolving 
dark energy models (i) by comparing the results of the 
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Union, the Constitution data sets, and the same Union 
members included in the Constitution sample and (ii) 
by analyzing the ACDM based mock data sets mimick- 
ing the Constitution and Union2 samples. 

According to our double-transition model, we found 
that the Union sample and the same Union members 
in the Constitution sample, both of which are based on 
the same light-curve fit parameters, show different be- 
haviors of dark energy equation of state due to different 
calibration processes applied to the distance moduli of 
the Union members. From the ACDM mock data anal- 
ysis, we also found that inhomogeneous distribution of 
SNIa data points in the redshift space and inhomoge- 
neous distance modulus errors cause biases in the dark 
energy equation of state parameters. We have demon- 
strated that the peak of the distribution of arithmetic 
means (estimated from thousands of mock data sets) can 



where w(z) is the dark energy equation of state param- 
eter at redshift z; z = ao/a — 1 with a(t) the cosmic 
expansion scale factor; H = a/a is the Hubble param- 
eter; a dot represents a derivative with respect to the 
cosmic time t; fi r , fi m , and fiDE are radiation, mat- 
ter, and dark energy density parameters, respectively; 
fix = —K/{a 2 H 2 ); subscript indicates present epoch. 
We set the present CMB temperature Tq = 2.725 K and 
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where wt (i = 0,1,2) is the dark energy equation of 
state parameter in each redshift bin; z trl and z tr2 rep- 
resent the locations of transition. The general formula 
with arbitrary number of bins can be found in Rcfs. 
[28l |33| . We assume a fiat background fi k = 0, thus 
fiDE = 1 — fir — fim- As we will consider models with 
fixed transition redshifts throughout this paper, we have 
five free parameters for our double-transition model. We 
denote the parameters as 6 = (h,£l m o,Wo,wx,W2); h 
is a normalized present day Hubble parameter Hq = 
lOO/i km s" 1 Mpc" 1 . 



be used as the unbiased estimator for the average bias in 
the dark energy equation of state parameters. 

The paper is organized as follows. We introduce our 
dark energy models in Sec. |TT] and describe our analysis 
method and the observational data used in this study in 
Sec. IIII1 Our results are presented in Sees. IIVI and IVl 
Section IVII is a discussion of our analysis with conclu- 
sions. 



II. DARK ENERGY MODEL 

We consider a wCDM model dominated by dark energy 
fluid with equation of state parameter w and cold dark 
matter. For generally varying dark energy equation of 
state the Fricdmann equation is given by 

(1) 

(2) 

I 

the number of massless neutrino species N y = 3.04. In 
this work we consider dark energy models with piece- 
wise constant dark energy equation of state, in which w 
parameter is constant for given redshift bins with a sud- 
den jump at transitions. As a simple case, we consider 
a double-transition model (with three redshift bins) . By 
imposing conditions that a and a are continuous at the 
transition epochs, we rewrite Eq. (J2J as 
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III. METHOD AND DATA 



A. Markov Chain Monte Carlo method 



To obtain likelihood distributions for parameters 6, wc 
use the Markov Chain Monte Carlo (MCMC) method. 
The MCMC chain elements are generated to ran- 
domly explore the whole parameter space based on the 
Metropolis-Hastings algorithm [46I . |47I ] . For the proba- 
bility density function P(6\D) which is needed to make 
decisions for accepting/rejecting a randomly chosen chain 
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element, we use 

Dfa\rt\ f 2 \^ [Xi(9) - X ohs l ] 2 
P(0|D) oc exp I -— J , x =2^ 2 » 

V Z / i= i °obs,i 

(4) 

where D denotes data, Xi{9) represents the model pre- 
diction for the iih observed data point X b s ,i with mea- 
surement error <7 bs,i, and N is the total number of data 
points. We use SNIa, BAO, and CMB data to constrain 
our model. In this case the final \ 2 is the sum of indi- 
vidual x 2 's: X 2 = Xsn + xIao + Xcmb- 

In applying the MCMC method to our dark energy 
models, we impose two priors. First, we impose a prior 
that the dark energy at sufficiently early epoch (zi = 10 8 ) 
should be lower than the maximum amount allowed by 
the big bang nucleosynthesis calculation, ^dej < 0.045 
[48j . If the number and the quality of SNIa data points 
within a redshift bin are not sufficient, the MCMC chain 
of the dark energy equation of state parameter of the 
bin usually does not converge, gradually decreasing to 
smaller and smaller values, sometimes reaching —100 in 
an extreme case. Therefore, as the second prior we set 
lower bounds to all the dark energy equation of state 
parameters, Wi > —5, to make sure the convergence of 
the MCMC chain. 

In applying the MCMC method to the observational 
data, we generate one million chain elements after a burn- 
ing period during which the first two thousand chain el- 
ements are discarded. For a ACDM based mock data set 
(see Sec. MI Cl below), we generate 0.2 million chain ele- 
ments. To test the convergence of MCMC chains we use 
a simple diagnostic: the means calculated from the first 
(after burning process) and the last 10% of the chain be- 
comes approximately equal to each other if the chain has 
converged (see Appendix B of Ref. [49| for detailed de- 
scription). We found that the 0.2 million chain elements 
are sufficient for the convergence of our MCMC chains. 
For each dark energy equation of state parameter in our 
models, we obtain its one-dimensional likelihood distri- 
bution by marginalizing over all other parameters, and 
compute the one-dimensional likelihood peak value or the 
arithmetic mean value as our average estimate, together 
with lcr (68.3%) and 2cr (95.4%) confidence limits. 

It should be noted that in our simple piecewise con- 
stant parametrization of dark energy equation of state 
we do not apply the decorrelating technique to obtain 
the independent estimate of w^s. Instead, by comparing 
the Wi's measured from the observational data with the 
likelihood confidence regions that the ACDM mock data 
analysis allows, we obtain a desirable statistical criterion 
for whether the measured deviation can be interpreted 
as dark energy evolution or not. 

B. Observational data 



BAO A, and CMB R parameters. 

For the SNIa data, we use the Union sample (307 mem- 
bers) 0, the Constitution sample (397 members) [||, and 
the Union2 sample (557 members) For Union and 
Union2 data analysis, we use the covariance matrix be- 
tween distance modulus errors without systcmatics. The 
Constitution data set consists of the same SNIa mem- 
bers in the Union data set and 90 CfA3 SNIa added at 
low redshift (z < 0.1). We will analyze four data sets: 
the Constitution set, the Union set, pure Union sample 
present in the Constitution set (hereafter Constitution- 
U), and Union2 set. The SNIa distance modulus in the 
Constitution-U set is different from the one in the original 
Union set due to a new calibration used in the Constitu- 
tion data set. 

Authors of Ref. Q produced the Constitution data set 
by taking directly the SALT light-curve ntparameters of 
the Union set (listed in Table 11 of Ref. [7|) and adding 
them to the SALT output of CfA3 sample. They found 
the best-fit cosmology from the combined SALT output 
parameters (rest- frame peak magnitude mj, stretch s, 
color c) by using the empirical relation for the distance 
modulus 

\i = m B - Mb + a(s - 1) - /3c, (5) 

obtaining empirical coefficients related to stretch correc- 
tion (a = 1.34+^|), color correction (ft = 2.59±g;J|), 
and absolute magnitude (Mb), which are different from 
those used in the Union data analysis (a = 1.24 ± 0.10, 
P = 2.28 ±0.11; 0). Furthermore, in the Constitution 
data set, a host galaxy peculiar velocity of 400 km s _1 has 
been assumed (rather than 300 km s _1 used in the Union 
data set 0), and the distance modulus uncertainties of 
Union SNIa members have been modified to reproduce 
the same amount of uncertainty in the measurement of 
w (based on the constant w dark energy model) as in Ref. 
0. Besides, an uncertainty of 0.138 mag has been added 
in quadrature to the uncertainties of SALT light-curve fit 
for each CfA3 SNIa. Since different empirical coefficients 
and different level of additional uncertainties were used to 
calculate the distance modulus of each SNIa, the Union 
and the Constitution sets have experienced different cal- 
ibration process (see Ref. Q for a complete description) . 
Figurc[T]compares distance moduli of Constitution-U and 
Union SNIa samples, where for each SNIa member the 
mean-subtracted difference between two distance moduli 
divided by the measurement uncertainty is shown. Al- 
though distance modulus values from two samples do not 
show special trend, the plot shows a scatter in both pos- 
itive and negative directions, implying that SNIa mem- 
bers of Constitution-U and Union samples experienced 
the different calibration process. 

For a given background world model the distance mod- 
ulus at redshift z is given by 



To constrain our dark energy model we jointly use the 
three kinds of observational data such as the SNIa data, 
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where is the redshift of recombination. From the 
WMAP 5-year data we have JH] 



z* = 1090, R ohs = 1-710 ± 0.019. 



The x 2 of the R parameter is 
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FIG. 1: A plot showing the difference between Constitution- 
U (jtic-u) and Union (/xu) SNIa distance modulus relative to 
the measurement error. Since each data set has an arbitrary 
zero point in the distance modulus, actually shown are mean- 
subtracted differences between two distance moduli divided 
by the distance modulus uncertainty of the Union SNIa (<7u)- 



with the comoving distance at redshift z given as 
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The SNIa data provide a distance modulus /i bs,i with 
a measurement error er bs,i for each SNIa at redshift Zi. 
We obtain x 2 for the SNIa data as 
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where iVgN is the number of SNIa in the data set, and 
^i{0) is the distance modulus at redshift Zi predicted by 
the theoretical world model to be constrained. 
The BAO A parameter is defined as [42[ 
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From the Sloan Dig ital Sky Survey luminous red galaxy 
sample we have [42| 
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where n s is the spectral index of scalar- type perturbation. 
We take n s = 0.960 based on the Wilkinson Microwave 
Anisotropy Probe (WMAP) 5-year results [H]. The \ 2 
of BAO data is 
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with a ohsA = 0.017 [42J. 

We use CMB shift parameter R introduced as [3, 145 
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C. Mock data 

In our estimation of Wi at each redshift bin, the devi- 
ation from ACDM value can be interpreted as the evo- 
lution of dark energy. However, even the observational 
data obtained in a perfect ACDM world may give a spu- 
rious deviation from w = — 1 due to inhomogencous dis- 
tribution of data points in the redshift space and the 
measurement errors, which makes the detection of dark 
energy evolution more complicated. To include such a 
spurious statistical effect in our analysis, we have gener- 
ated thousands of ACDM-motivated mock data sets and 
analyzed them in a similar way as the real observational 
data is analyzed. 

For SNIa mock data, we make the five thousand Con- 
stitution and Union2 mock data sets predicted in the 
ACDM model by using rcdshifts (zi's) and distance mod- 
ulus measurement errors (<J bs,i's) of SNIa members in 
the Constitution and the Union2 samples. For ACDM 
background world model we use cosmological parameters 
that are consistent with the WMAP 5-year observation, 
i.e., h = 0.705, Q, mQ = 0.274, fl K0 = 0, and w = -1 @. 
For a given SNIa at redshift z%, the mock data point is 
generated by ^ moc k,i = Macdivu + &i, where /iacdm,; is 
a distance modulus predicted in the ACDM model and 
ej is a random number drawn from a Gaussian normal 
distribution with variance <J^ hs ., which mimics the mea- 
surement noise. In this manner, we can get ACDM based 
SNIa mock data sets with the same redshift distribu- 
tion and distance modulus precision of Constitution and 
Union2 samples. The same method is applied to generat- 
ing mock BAO A and CMB R parameters. The similar 
technique of gen erating SNIa mock data sets has been 
used in Refs. [13, HJ. 

There are two good points of using ACDM mock data 
sets. First, we can include the statistical effect induced 
by different random realizations of measurement errors, 
which usually appears as additional variance in our es- 
timation of dark energy equation of state parameters. 
Thus, the confidence regions of Wi's obtained from the 
ACDM mock data sets become bigger than those from 
the single data set (see Figs. [3] and S] below). By com- 
paring results from the observational data set with those 
from ACDM mock data sets, we can avoid a false signa- 
ture of dark energy evolution due to statistically peculiar 
measurement noise in the observational data. Analyzing 
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the single data set cannot account for such a statistical 
effect. Second, we can estimate any bias in the estimate 
of Wi relative to — 1, induced by the sparse number of 
data points and measurement errors within the redshift 
bin considered. We investigate this issue in detail in Sec. 

El 



IV. CONSTRAINTS ON DARK ENERGY 
EQUATION OF STATE IN 
DOUBLE-TRANSITION MODELS 

In this section, we present the results of our dark en- 
ergy model with double transitions. Figure [2] shows Wo, 
wi, and W2 measured from Constitution, Constitution- 
U, Union, and Union2 SNIa samples together with BAO 
and CMB data, for various transition rcdshifts. We fix 
the width between the two transitions as Ztr2 — ^tri = 0.2, 
and vary zm from 0.1 to 0.5 with steps of 0.1. From the 
marginalized one-dimensional likelihood distributions of 
u>o, w±, and W2 parameters, we have estimated peak lo- 
cations (denoted as P) and la (68.3%), 2a (95.4%) confi- 
dence limits (points with error bars) . In the Constitution 
and Constitution-U sets wq is larger than —1 and w\ is 
smaller than —1 for z trl < 0.4. However, such trends 
are not seen in the Union and Union2 sets. Values of u>2 
tend to be larger than —1 in all four data sets. Since the 
Constitution-U sample is composed of exactly the same 
SNIa members of the Union data set, we can conclude 
that such different behaviors of Wq and wi arise not due 
to the new CfA3 data added in the Constitution set but 
due to the different calibration applied to each data set. 

Now, we consider a case of z tr i = 0.2 and z tT 2 = 0.4 
where the deviations in Constitution and Constitution- 
U sets are most significant. The result is shown in the 
top row of Fig. [3] In the Constitution and Constitution- 
U sets the la confidence limits of all w^s deviate from 
w = — 1. However, in the cases of Union and Union2, 
Wq and Wi are consistent with ACDM whereas la limits 
of u>2 deviate from w = — 1. In order to estimate the 
statistical significance of the trends, we have also ana- 
lyzed ACDM based Constitution mock data sets, and es- 
timated 68.3% and 95.4% confidence limits on Wi's based 
on the total sum of five thousand MCMC chains; these 
are shown as grey regions in Fig. [3] (top panels). The 
likelihood confidence regions for the ACDM mock data 
sets are wider than those determined from the marginal- 
ized one-dimensional likelihood distributions for the Con- 
stitution data set (red error bars). This is because the 
mock data sets include the additional statistical variance 
caused by different realizations of measurement errors. 

As a second measure of average, we compute the arith- 
metic mean (denoted as M) of the MCMC chain elements 
for each cosmological parameter. In fact, it is known that 
in the MCMC method the arithmetic mean of MCMC 
chain is appropriate as the best estimate for the param- 
eter and is a more robust quantity than the peak (or 
mode) of the marginalized distribution (52|. The bot- 



tom row of Fig. [3] shows distributions of five thousand 
M values of wq, W\, and u>2, each estimated from the 
MCMC chain of individual mock data set. The smooth 
distributions (shown as dark red curves) have been gen- 
erated by applying the Savitzky-Golay smoothing filter 
[53l ] to the histograms with the bin size Aw = 0.02. The 
distribution for wi is wider than those for wo and u>2, 
implying that the equation of state parameter w\ is less 
tightly constrained in our double-transition model with 
z trl =0.2; this is due to different number of data in each 
bin. For a comparison, M values for the observed SNIa 
data sets are indicated by vertical lines. In Constitution 
and Constitution-U sets all M values of w;'s are outside 
la regions of the distributions. In the Union set, the M 
values of wq and w\ are consistent with the ACDM mock 
data, whereas that of u>2 is outside the la region. Al- 
though we detect more than la deviations of all Wi's in 
the Constitution and the Constitution-U sets, we would 
like to emphasize that all the deviations are still within 
2a regions, consistent with ACDM model. The similar- 
ity between Constitution and Constitution-U sets and the 
difference between Constitution-U and Union are impor- 
tant because the Constitution-U and the Union sets are 
based on the same SALT light-curve fit parameters of the 
Union SNIa sample. The different behaviors of wq and 
Wi between the Constitution-U and the Union sets are 
purely due to different calibration processes applied in 
producing the SNIa distance modulus (see Sec. IIII 51) . 

Table U compares equation of state parameters and 
matter density parameter in the double-transition dark 
energy model with z t ri = 0.2 and Za-2 = 0.4 for vari- 
ous SNIa data sets and average estimation methods. For 
ACDM based Constitution mock data sets, first we ob- 
tain distributions of P or M values and then present 
the peak of P distribution (PP), the peak of M distri- 
bution (PM), and the mean of M distribution (MM), 
together with 68.3% confidence limits. The location of 
peak and the confidence limits are determined based 
on the smooth distribution generated by applying the 
Savitzky-Golay smoothing filter to the histograms of P 
or M values. Comparison of the parameters estimated 
by three different average methods (PP, PM, and MM) 
with the fiducial value w = — 1 suggests that wi is biased 
to the positive direction by the amount of about 0.07- 
0.09. The estimated bias is uncertain roughly within 
a/^N ~ 0.45/V5000 = 0.006. The w is very weakly 
biased to the negative direction. For u>2 parameter, the 
presence of bias is not clear at this point because the three 
average estimates PP, PM, and MM indicate positive, 
zero-consistent, and negative biases relative to w = — 1, 
respectively. 

We notice all the observational data sets show a trend 
that P value is larger than M value for W2 parameter, 
u>2 (P) > W2 (M) unlike other parameters. The similar 
feature is seen in the results of Constitution mock data 
analysis (PP and PM values in Table H]). As shown in 
bottom row of Fig. [31 in the case of Constitution mock 
data sets M distributions for Wq and wi look symmetric 
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FIG. 2: Estimated Wo, ttfi, and W2 in the double-transition dark energy models. The transition redshift z tI i is varied from 0.1 
to 0.5 by 0.1 increase in redshift while the width between the two transitions is fixed as zti2 — Ztri = 0.2. Red, green, blue, and 
brown colors represent Constitution, Constitution-U, Union, and Union2, respectively. The error bars with dark (light) color 
indicate 68.3% (94.5%) confidence limits. 




FIG. 3: Top panels: Estimated wo, Wi, and u>2 in the double-transition dark energy model with z tr i = 0.2 and jz t r2 = 0.4, 
taken from Fig. [2] with the same colored code. Data points with error bars are based on the actual SNIa data sets. In all 
panels, shaded regions indicating 68.3% (dark grey) and 95.4% (light grey) confidence limits have been obtained from the five 
thousand ACDM based Constitution mock data sets. Bottom panels: Distributions of arithmetic mean (M) values of wo, Wi, 
and W2 obtained from the Constitution mock data sets (histograms fitted by dark red smooth curves, see text). The three bins 
in the Constitution set have 147, 49, and 201 SNIa data points, respectively. Red (solid), green (dotted), blue (dashed), and 
brown (dot-dashed) vertical lines denote locations of M values for Constitution, Constitution-U, Union, and Union2 data sets, 
respectively. 
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TABLE I: Estimated u>o, wi, W2, and f2 m o in the double-transition dark energy model for z tr i = 0.2 and z tr 2 = 0.4. For 
comparison, results from the five thousand ACDM based Constitution mock data sets are presented in the bottom. Symbols 
in the parenthesis represent our methods of average estimation, and the errors indicate 68.3% confidence limits (see text). 
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u ' z ' ' -0.017 

n 970+0.014 
U.z!C8_ 017 
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while that for W2 has a non-symmetric shape with a tail in 
the direction to which W2 decreases. We expect from the 
PP value together with ±ler errors listed in Table Q] that 
the P distribution for wi is also asymmetric and has a 
tail in the same direction and a peak location larger than 
that of the M distribution. In the increasing direction of 
W2 the parameter is constrained by the big-bang nucle- 
osynthesis calculation and the CMB shift parameter. On 
the other hand, W2 is free to randomly walk in the op- 
posite direction without a strong prior like u>2 > —5 (see 
Sec. IIII[) since the subdominant dark energy component 
at high redshift has more freedom to take any value of 
equation of state parameter. That is the reason why the 
likelihood distribution of W2 parameter has a tendency of 
having the asymmetric shape with a tail to the smaller 
u>2 direction. Thus, it is naturally expected that the lo- 
cation of peak (P) in the one-dimensional marginalized 
distribution is larger than the arithmetic mean (M) of 
MCMC chain elements. 

We made similar analysis based on Union2 mock data 
sets, and the results are shown in Fig. @] and Table [TT1 
Here, we choose the transition redshifts z tr i = 0.1 and 
zti-2 = 0.3, for which the deviations of Wi's from w = — 1 
are the largest in the case of Union2 data set (see Fig. [2]) . 
In the case of Union2 set, la confidence limits of w\ and 
W2 deviate from w = — 1 whereas wq is consistent with 
ACDM. However, all the values of w^s from the Union2, 
Union, Constitution, and Constitution-U sets are within 
2(j limits thus are consistent with ACDM model (top row 
of Fig. [4j. We also see that compared with the M dis- 
tributions for ACDM based Union2 mock data sets all 
M values of w^s are consistent with ACDM model and 
deviations are not statistically significant. The M distri- 
butions of W\ and W2 (bottom row of Fig. [4]) are narrower 
than those for Constitution mock data sets (bottom row 
of Fig. [3]), reflecting that the dark energy equation of 
state parameters are relatively well constrained as we 



have more SNIa data points in the redshift bins. 

We observe the same trend that P value is larger than 
M value for w 2 parameter, w 2 (P) > W2 (M). But the 
difference between P and M values is smaller due to the 
increase of the width of the final redshift bin (z > 0.3) 
and of SNIa data points within the bin. As shown in 
Table UT1 the Union2 mock data analysis implies that es- 
timated wo and Wi are unbiased by three average meth- 
ods (PP, PM, and MM); they are almost the same as the 
fiducial value w = — 1. The PP, PM, and MM average 
methods for W2 parameter give different values, indicat- 
ing the likelihood distribution for this parameter is also 
not symmetric. However, the asymmetry is weaker than 
the case of Constitution data set. 



V. CHARACTERS OF CURRENT SNIA DATA 
BASED ON MOCK DATA ANALYSIS 

In the double-transition model the Constitution mock 
data sets shows some statistical deviation of equation of 
state parameters from the fiducial ACDM model (Table 
U). For PM average, the value Wi = —0.92 indicates a 
bias at intermediate redshift interval although it is small 
compared with the size of la error. In the Constitution 
sample, the SNIa data points are irregularly distributed 
between z = 0.015 and z = 1.551 with non-uniform dis- 
tance modulus errors, and the data is sparse in the inter- 
mediate redshift z = 0.1-0.3. In this section, we investi- 
gate a bias that may be induced by the characters of the 
current SNIa data sets. The estimation of the bias is im- 
portant because it can be misinterpreted as the spurious 
evolution of dark energy. 

We have generated ACDM based mock data sets in 
three different ways. First, we made SNIa data with 
uniform distribution in the redshift interval and uniform 
distance modulus errors (Case 1). We arrange SNIa 
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FIG. 4: The same as Fig. [3] but for z tI i =0.1 and Ztr2 = 0.3 and based on Union2 mock data sets. The three bins in Union2 
set have 165, 117, and 275 SNIa data points, respectively. In the top panels, data points with error bars are based on the actual 
SNIa data sets while the shaded regions (in all panels) are obtained from the five thousand ACDM based Union2 mock data 
sets. 



TABLE II: The same as Table U but for z tI i = 0.1 and Zt r 2 = 0.3 and based on Union2 mock data sets. 



.020 
.017 
+0.019 



Data Average wo wi W2 £lm0 

Constitution (P) -0.53loJl -1.33±g;^ -0.80±^f 0.2851°,;°,^ 

Constitution-U (P) -0.77^1 -1.22l°^ -0-82+g^ 0.2801°, "~ 

Union (P) -0.95t°:g -1.10±S;g -0.74t°:JS 0.280l°. ol7 

Union2 (P) -0.80tg;l -1.26±g;l° -0.77±g-^g 0.278±g;gjg 

Constitution (M) — 0.56±g;|J -1.32±°;^ -0.861°,;?!? 0.2861°,;°,^ 

Constitution-U (M) -0.70lg g — 1.24±g-|g -0.88±g"^ 0.2831^14 

Union (M) -0.98±° o f 8 -1.14±g:» -0.79t°;? 2 5 0.2821°;°" 

Union2 (M) -OSOlgg -L27lgg -O-Sllgg 0.280±g;g^ 

Union2 mock (PP) -1.001°,;^ -1.001°,;^ — 0.95±g;|| 0.276±g:g^ 

Union2 mock (PM) -l.OOl ,;^ -1.001°,;^ -l.Ol+g;^ 0.2781°;°^ 

Union2 mock (MM) -l.Oll^i -1-OOlg;^ -1.06±g;?? 0.278±g;gi| 



data points (ATsn = 397) uniformly on the rcdshift range 
spanned by the Constitution sample and assign the uni- 
form error to all members. We take the harmonic average 
as the uniform error given by 

a = J N NSN - = 0.177, (15) 

in magnitude scale. The second-type mock data sets have 
the uniform redshift-distribution but with distance mod- 
ulus errors of the Constitution sample (Case 2). We sort 
the Constitution SNIa data points in increasing order in 



rcdshift, then assign distance modulus errors to mock 
SNIa members in the same order. The third-type mock 
data sets have the same redshift-distribution as the Con- 
stitution data set but with uniform distance modulus er- 
rors (Case 3). For comparison, we include the Constitu- 
tion mock data sets as Case 4; these are already presented 
in Table |U We have analyzed these mock data sets (in- 
cluding mock BAO and CMB parameters) in the same 
way as the Constitution mock data sets have been an- 
alyzed for double-transition model with z t ri = 0.2 and 
z tT 2 = 0.4 (Sec. |TVj). The results are shown in Table [TTTl 
which compares dark energy equation of state parameters 
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obtained by using different average estimation methods 
(PP, PM, and MM). 

For all Cases 1-4, we see that wq and W\ almost 
do not depend on the choice of average methods, that 
is, w l (PP) ~ w t (PM) ~ Wi (MM), while w 2 pa- 
rameters show a trend shown in the previous section, 
w 2 (PP) > w 2 (PM) > w 2 (MM). The likelihood distri- 
butions of wq and w\ are symmetric whereas that of w 2 
is not symmetric for all cases of mock data sets. This is 
verified in Fig. [5l which compares distributions of M val- 
ues derived from the four cases of ACDM mock data sets. 
Expecting that analysis of Case 1 (uniform distribution 
and error) restores the fiducial value w = — 1 correctly, 
we found that the PM acts as the unbiased (or less bi- 
ased) estimator for the average of parameter because PM 
values of wq , w\ , and W2 are consistent with the ACDM 
model. Even for the non-symmetric shape of likelihood 
distribution as in the case of W2 parameter, the PM aver- 
age is unbiased, and thus provides the precise estimation 
of bias in the parameter. 

For Constitution mock data sets (Case 4), the bias in 
the positive direction for wi parameter is observed. Be- 
cause the mock data sets of Case 3 (uniform error) induce 
the similar amount of bias in the same direction while the 
data sets of Case 2 (uniform distribution) does not affect 
wi much, we conclude that the irregular distribution of 
Constitution SNIa data points is more concerned with 
the bias than the non-uniform distance modulus errors. 
However, the bias estimated based on the ACDM mock 
data sets is small compared with the precision the cur- 
rent SNIa data (including BAO and CMB parameters) 
allows. For Constitution mock data (Case 4), it amounts 
to 1 + Wi ~ 0.18cr at most for PM average. Unlike the 
bias in the w\ parameter, W2 is not much biased (based 
on PM average). 

For the Union2 mock data sets, we have not seen any 
significant deviation of equation of state parameters (see 
Table HI)) . Because of the relatively denser SNIa data 
points, the biases in Wi estimation are smaller than those 
of Constitution case. However, the likelihood distribu- 
tion of w 2 still appears asymmetric (as shown in Fig. 
[4]), and thus the estimated value of W2 depends on the 
average methods. 



VI. DISCUSSION 

In this study, we have investigated the evolution of 
dark energy equation of state parameter w(z) using the 
sudden-jump approximation of w for some chosen red- 
shift intervals with double transitions. We used four SNIa 
data (Constitution, Constitution-U, Union, and Union2) 
together with BAO A parameter and CMB R parameter, 
and used the MCMC method with Metropolis-Hastings 
algorithm to obtain the likelihood in the parameter space. 
In order to test the statistical significance of, or to esti- 
mate the bias in, the equation of state parameters we 
investigated the case of ACDM mock data sets. Since 



the analysis of thousands of ACDM mock data sets can 
include all statistical effects due to the characters of data 
sets such as the irregular distribution in the redshift inter- 
val and the different random realization of measurement 
errors, our method of using the ACDM mock data is com- 
plementary to the uncorrclatcd estimation of dark energy 
equation of state parameters based on the decorrelating 
technique [H, [H, 111]. 

In the double-transition dark energy models, we show 
that there are some deviations from ACDM in the low 
(z < 0.2) and the middle (0.2 < z < 0.4) redshift re- 
gions especially in the Constitution and Constitution-U 
sets. Deviations of Wi's from ACDM in the low and the 
middle redshift are between 1 and 2a confidence levels, 
see the top row of Fig. [31 The deviations are stronger 
when we compare each arithmetic mean (M) value of Wi 
with the ACDM mock data; the M values of Constitution 
and Constitution-U sets are situated near 2a boundary 
of the ACDM mock data distribution, see the bottom 
row of Fig. [31 whereas there are no such deviations in the 
Union set. Thus, we conclude that our analysis of the 
double-transition model indicates statistically noticeable 
deviation of the Constitution data from ACDM model in 
the low redshift. Furthermore, consistent result between 
Constitution set and Constitution-U set suggests that the 
detected deviation in the Constitution data is not due to 
additional data in the Constitution data, but mainly due 
to different calibration of the same Union data. Results 
of Constitution-U compared with the Constitution and 
the Union data in Figs. [2] and [3J apparently show that 
the new 90 CfA3 SNIa data added in the Constitution 
data do not have role in causing the difference between 
the Constitution and the Union sets. 

Our conclusion is still valid when we measure dark en- 
ergy equation of state parameters by choosing different 
transition redshifts, as summarized in Fig. [5] (top three 
rows), where the width of the second redshift bin has 
been extended by the amount of Az = 0.1 to left, right, 
and both directions. Although extending the redshift in- 
terval of the middle redshift bin makes the error bar of 
wi value smaller as expected, the observed trend for wq 
and wi estimates of Constitution, Constitution-U, and 
Union samples docs not change much in comparison with 
the result shown in Fig. [3J 

In the bottom row of Fig. [S] we present our measure- 
ment of equation of state parameters based on the triple- 
transition model with the first three redshift bins equally 
spaced between z = and 0.4. Even in this case our con- 
clusion maintains. However, increasing number of red- 
shift bins causes error bars to be drastically larger and 
it is more and more difficult to make the MCMC chains 
converge without strong priors. This is why we have cho- 
sen a simple dark energy model with double transitions. 

We emphasize that the Union and the Constitution- 
U sets are composed of the same SNIa members and 
arc based on the same light-curve fit parameters. In a 
dark energy model with constant equation of state, the w 
parameter (P value) estimated from SNIa data together 
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FIG. 5: Distributions of M values for four different types of ACDM mock data sets. The smooth distributions have been made 
based on the histograms with the bin size Aw — 0.001. Black (solid), blue (dashed), yellow (dot-dashed), and red (dotted) 
curves represent Case 1 (uniform distribution & error), Case 2 (uniform distribution), Case 3 (uniform error), and Case 4 
(Constitution mock), respectively. Vertical dashed lines indicate the location of w — — 1. 
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FIG. 6: Top three rows: dark energy equation of state parameters (too, Wi, W2) measured from Constitution, Constitution-U, 
Union, and Union2 samples for different choices of transition redshifts, ztri and ztr2- In order from top to bottom we choose (a) 
ztri =0.1 and z tr 2 = 0.4, (b) z tI i = 0.2 and z tI 2 = 0.5, and (c) z t ri =0.1 and z tI 2 = 0.5. Bottom row: (d) dark energy equation 
of state parameters in the triple-transition dark energy model. The first three bins are equally spaced in the redshift interval 
2 = 0-0.4. Each number listed in each redshift bin indicates the number of SNIa included in the bin. Note that although the 
size of error bars generally gets smaller as the corresponding redshift bin includes more SNIa data points, it is also affected by 
the quality of SNIa data points within the bin. 
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TABLE III: Dark energy equation of state parameters of double-transition model with ztti = 0.2 and Ztr2 = 0.4 estimated from 
the five thousand ACDM based SNIa mock data sets with four different types (Cases 1-4; see text). Values are listed based on 
three average estimation methods (PP, PM, and MM). The errors indicate 68.3% confidence levels. 



Data 



!('() 



(PP) Wl (PP) W 2 (PP) 



Case 1 
Case 2 
Case 3 
Case 4 



uniform distribution and error —0.99 

uniform distribution —0.97 

uniform error —1.01 

Constitution mock —1.02 



+0.27 
0.27 
+0.29 
0.30 
+0.18 
0.19 
+0.17 
0.17 



-1 00 +u ' 4S 

1 - uu -0.47 
.1 no+0.46 
1 - UO -0.47 
+0.46 
0.50 
+0.43 
0.47 



-0 



+0.22 
-0.25 

0.95_ a22 



-0.97 
-0.93 



-0.25 
-,+0.27 
-0.34 

0.92-;:;- 



-0.92: 



Data 



w (PM) Wl (PM) w 2 (PM) 



Case 1 
Case 2 
Case 3 
Case 4 



uniform distribution and error 
uniform distribution 
uniform error 
Constitution mock 



-0 



l.UU_ js 
+0.30 
-0.31 

-1 oi+ 017 

l-Ul_o.l 9 

-i n3 +017 



+0.48 
0.47 
+0.45 
-0.49 

u ' yd -0.49 
+0.46 
0.45 



-1.00 
-1.01 



-0.92 



-1.01 
-1.03 
-1.02 
-1.02 



+0.24 
0.29 
+0.25 
0.27 
+0.24 
0.27 
+0.28 
0.47 



Data 



wo (MM) iui (MM) w 2 (MM) 



Case 1 
Case 2 
Case 3 
Case 4 



uniform distribution and error 
uniform distribution 
uniform error 
Constitution mock 



i nn +0 - 28 



-0.98 
-1.03 
-1.03 



+0.30 
0.31 
+ 11.1!) 

0.17 
+0.17 
0.17 



_1 (11+0.49 

-1 04+ - 48 
-0 94+ 47 

U-»*_0.48 

-0 91+ - 45 



-1.10 
-1.08 

-1 

-1.23 



+0.33 
0.20 
+0.30 
0.22 
-^■7+0.39 



0.12 
+0.07 
0.26 



with BAO A parameter is l+w = -O.Oll+g^so for Union 
sample and l+w = — 0.010_Q'ggg for Constitution-U sam- 
ple, both of which are nearly identical to each other and 
arc similar to the values obtained by Refs. @, _}. Al- 
though the two data sets give nearly identical results in 
the simple constant w dark energy model, they show dif- 
ferent behaviors of dark energy equation of state in our 
evolving dark energy model with double transitions, due 
to different calibrations applied during the production of 
SNIa distance modulus. Our analysis suggests that it is 
generally important to distinguish the real w evolution 
from the calibration artefact in the SNIa data analysis. 

However, we cannot conclude that such deviations in 
the Constitution data are the strong hint for evolution of 
dark energy because all the equation of state parameters 
Wi's measured from all the observational data sets are 
consistent with ACDM model within 2a confidence lim- 
its. Our results are also consistent with previous studies 
searching for dark energy evolution with SNIa data based 
on the redshift- binning method 0, [3^-[__| . Furthermore, 
the difference between Union and Constitution-U due to 
different calibrations is also not significant at the level of 
current precision of SNIa data. Such an effect caused by 
different calibrations should be considered to be impor- 
tant in the forthcoming missions of SNIa survey. 

From the analysis of ACDM-motivatcd mock data sets 
with four different types (Cases 1-4), we have tried to 
estimate bias in the equation of state parameters due to 
the character of data set used in search for dark energy 
evolution. Comparing average values for Wi's obtained 
by three different methods such as PP (peak of peaks), 



PM (peak of means), and MM (mean of means) while 
knowing the answer (w = —1), we conclude that the PM 
is an unbiased estimator of average bias. For the sym- 
metric likelihood distribution, all average methods give 
values consistent with each other (e.g., see Fig. [3] and 
wq and wi for the Constitution data sets in Tabic IIIII) . 
However, for the non-symmetric likelihood distribution 
as in the case of «_ parameter, the PP and MM aver- 
ages are biased to the opposite directions with respect to 
the PM average which restores the fiducial value w = — 1 
correctly. The Constitution sample turns out to induce 
a bias in the wi parameter by the amount of about 0.2a 
whereas the Union2 sample does not show any noticeable 
bias in the estimation of parameters (see Table [H]) . In 
order to precisely estimate any bias due to the character 
of data set based on the mock data analysis, we suggest 
to use the PM quantity, the peak (or mode) of the dis- 
tribution of arithmetic means, as the unbiased estimator 
for the average bias. 



Acknowledgments 

We wish to thank Dr. Eric V. Lindcr, Dr. Yong-Scon 
Song, and Professor Dal Ho Kim for useful comments 
and discussions. We also would like to thank the anony- 
mous referee for the constructive and helpful comments 
on our manuscript. This work was supported by the Ko- 
rea Research Foundation Grant funded by the Korean 
Government (KRF-2008-341-C00022). 



12 



A.G. Riess et al. As tron. J. 116, 1009 (1998) 

J. 517, 565 (1999) 



arXiv:astro-ph/9805201 



Perlmutter et aL Astrophys 
[arXiv:astro-ph/9812133] . 
A . G. Riess et al. Astrophys . 



J. 659, 98 (2007) 
Astier et al. Astron. Astrophys. 447, 31 (2006) 



arXiv:astro-ph/0611572 



arXiv:astro-ph/0510447 



T.M. Davis et al. As trophys. J. 666, 716 (2007) 
[arXiv:astro-ph/0701510] . 

H. La mpeitl et al. Mon. Not. R. Astron. Soc. 401, 2331 

(2009) |arXiv:0910.2T93l [astro-ph.CO]] . 

Supernova Cosmology Project Collaboration, M. Kowal- 

ski et al. Astrophys. J. 686, 749 (2008) |arXiv:0804.4142l 

[astro-ph]]. 

M. Hicken et al. Astrophys. J. 700, 1097 (2009) 
[arXiv:0901.4804l [astro-ph.CO]]. 

R. Amanullah, e t al. Astrophys. J. 716, 712 (2010) 
[arXiv: 1004. 17TT1 [astro-ph.CO]]. 

G. Efstathiou, Mon. Not. R. Astron. Soc. 301, 842 (1999) 
[arXiv:astro-ph/9904356] . 

D. Hut erer and M.S. Turner, P hys. Rev. D 64, 123527 
(2001) [arXiv:astro-ph/0012510l . 
J. Well er and A. Albrecht, Ph ys. Rev. 
(2001) [arXiv:astro-ph/0008314l . 



Lett. 86, 1939 



9 



D. Hut erer and G. Starkman, Ph ys. Rev. Lett. 90, 031301 
(2003) [arXiv:astro-ph/0207517'] . 

P.S. Corasaniti and E.J. Copeland, Phys. Rev. D 67, 

063521 (2003) [arXiv:astro-ph/0205544] . 

U. Alam, V. Sahni, T.D. Saini and A.A. Starobin- 

sky, Mon. Not. R. A stron. Soc. 354, 275 (2004) 

[arXiv:astro-ph/0311364] . 

R.A. Daly and S.G. Djorgovski, Astrophys. J. 597, 

(2003) [arXiv:astro-ph/0305197'] . 

A.G. Kim, E.V. Linder, R. Miquel and N. Mostek, 
Mon. Not. R. Astro n. Soc. 347, 909 (2004) 
|arXiv:astro-ph/0304509] . 

H.K. Jassal, J.S. Bagla and T. Padmanabhan, 
Mon. Not. R. Astro n. Soc. 356, Lll (2005) 
|arXiv:astro-ph/0404378] . 

R.A. D aly and S.G. Djor govski, Astrophys. J. 612, 652 

(2004) [arXiv:astro-ph /0403664] . 
C. Wetterich, Phys. Lett. B 594, 17 
arXiv:astro-ph/0403289] . 

Grav 



Y.G. Gong, Class. Q uant. 
|arXiv:astro-ph/0405446] . 

D. Huterer and A. Cooray, Phys. 
(2005) [arXiv:astro-ph /0404062] . 

E. V. Linder and D. Huterer, Phys 
(2005) [arXiv:astro-ph/0505330l . 
Y.G. Gong and Y.Z. Zhang, P hys. 
(2005) [arXiv:astro-ph /0502262] 
A. Shafieloo, U 
sky, Mon. Not, 



22, 2121 



(2004) 
(2005) 



Rev. D 71, 023506 



Rev. D 72, 043509 



Rev. D 72, 043518 



Alam, V. Sahni and A.A. Starobin- 
R. Astron. Soc. 366, 1081 (2006) 



|arXiv:astro-ph/0505329| . 
A. Shafieloo, Mon. Not.R . Astron. Soc. 380, 1573 (2007) 
|arXiv:astro-ph/0703034] . 

R. de Putter and E.V. Linder, Astropart. Phys. 29, 424 
(2008) [arXiv:0710.0"373l [astro-ph]]. 



[28 
[29 
[30 
[31 
[32 
[33 
[34 
[35 
[36 
[37 

[38 
[39 
[40 

[41 
[42 
[43 
[44 
[45 
[46 

I 47 

[48 

[49 

[50 

[51 
[52 
[53 



S. Sulli van, A. Cooray a nd D.E. Holz, JCAP 0709, 004 

(2007) larXiv:0706.3"730l [astro-ph]]. 

A. Hojj ati, L. Pogosian a nd G.B. Zhao, JCAP 1004, 007 

(2010) |arXiv:0912.4843l [astro-ph.CO]]. 

J.B. Jimenez, R. Lazkoz and A.L. Maroto, Phys. Rev. D 

80, 023004 (2009) |arXiv:0904.0433V 2 [astro-ph.CO]]. 

Y. Wang, Ph ys. Rev. D 80, 123525 (2009) 

|arXiv:0910.2492l [astro-ph.CO]] . 

A. Shafieloo, V. Sahni and A.A. Starobinsky, Phys. Rev. 
D 80, 101301 (2009) |arXiv:0903.5T4"Tl [astro-ph.CO]]. 
Q.G. Huang, M. Li, X.D. Li and S. Wang, Phys. Rev. D 
80, 083515 (2009) |arXiv:0905.0"797l [astro-ph.CO]] . 
R.G. Cai, Q. Su a nd H.B. Zhang, JCAP 1004, 012 (2010) 
|arXiv:1001.2"207l [astro-ph.CO]] . 

G.B. Zhao and X. Zhang, Phys. Rev. D 81, 043518 (2010) 
|arXiv:0908. 15681 [astro-ph.CO]] . 

Y.G. Gong, R.G. Cai, Y. Chen a nd Z.H. Zhu, JCAP 
1001, 019 (2010) |arXiv:0909.0596l [astro-ph.CO]] . 
P. Serra, A. Cooray, D.E. Holz, A. Melchiorri, S. Pan- 
dolfi and D. Sar kar, Phys. Rev. D 80, 121302 (2009) 

|arXiv:0908.3"T86l [astro-ph.CO]]. 

S. Wang, X.-D. Li, and M. Li. larXiv:1009.5837fr 3 [astro- 
ph.CO]. 

S. Qi, T. Lu and F.-Y. Wang, Mon. Not. R. Astron. Soc. 
398, 78 (2009) larXiv:0904.2832l [astro-ph.CO], 
M. Chevalli er and D. Polarski, In t. J. Mod. Phys. D 10, 
213 (2001) [arXiv:gr-qc/0009008] . 
E.V. Linder, Phys. R ev. Lett 
[arXiv:astro-ph/02085i2] . 
D.J. Eisenstein et al. [SPSS Collaboration 
633, 560 (2005) 



90, 091301 (2003) 



Astrophys. 



J. 633, 560 (2005) [arXiv:astro-ph/0501171 . 
E. Komatsu et al. [WMAP Collaboration], Astrophys. J. 
Suppl. 180, 330 (2009) |arXiv:0803.0"547l [astro-ph]]. 
J.R. Bond, G. Efstathiou and M. Tegmark, Mon. Not. R . 
Astron. Soc. 291, L33 (1997) [arXiv:astro-ph/9702100] . 
Y. Wang and P. Mukherj ee, Astrophys. J. 650, 1 (2006) 
[arXiv:astro-ph /060405T] . 

N. Metropolis, A.W. Rosenbluth, M.N. Rosenbluth, A.H. 
Teller and E. Teller, J. Chem. Phys. 21, 1087 (1953). 
W.K. Hastings, Biometrika 57, 97 (1970). 
R. Bean, S.H. Hansen, and A. Melchiorri, Phys. Rev. D 
64, 103508 (2001) [arXiv:astro-ph/0104162]; B.A . Bas- 
sett et al., JCAP 07. 007 f2008) [arXiv:0709.05261 . 
A. Abrahamse, A. Albrecht, M. Barnard, and B. Bozek, 
Phys. Rev. D 77, 103503 (2008) |arXiv:0712.2"879l [astro- 
ph]]. 

T. Holsclaw, U. Alam, B. Sanso, H. Lee, K. Heitmann, 
S. Habib and D. Higdon, Phys. Rev. D 82, 103502 (2010) 

larXiv: 1009.54431 [astro-ph. CO] . 

I. Antoniou and L. Perivolaropoulos, arXiv: 1007.4347 
[astro-ph.CO]. 

D.N. Spergel, et al., Astro phys. J. Suppl. 148, 175 (2003) 
[arXiv:astro-ph /0302209] . 

W.H. Press, S.A. Teukolsky, W.T. Vetterling, and B.P. 
Flannery, Numerical Recipes in Fortran. The art of sci- 
entific computing, 2nd eds. (Cambridge University Press, 
1992), p. 644. 



